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Figure 1: Overall system flow 
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Figure 2a: Step 110 Selection Procedure 
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Figure 2b: Corresponding pseudo code (When Applicable comments point to the steps in the 
flowchart) 

Sub CPredict() 

Dim PlanningHorizon, RouteHorizon, SafetyHorizon,LTHorizon As Integer 

Dim CostQ As Double 

DimMinCost As Double 

Dim LastCustomer As Integer 

Dim Minday, calcday As date 

Dim Count As Integer 

Dim Db As Database 

Dim Customers As Recordset 

Set Db = CurrentDb 

'400 Read information 

RouteHorizon = 7 

SafetyHorizon=2 

LTHorizon=7 

ReDim Cost(PlanningHorizon) 

Set Customers = Db.OpenRecordset("modelreg") 
PlanningHorizon = RouteHorizon+SafetyHorizon+LTHorizon 

'402 Select first customer 

Customers.MoveFirst 

Customers.MoveLast 

LastCustomer = Customers .RecordCount 

Customers.MoveFirst 

For Customer = 1 To LastCustomer 
'404 Initialization 
MinCost = 99999999 
calcday = Now 

'406 Calculate cost 

For Day = 1 To PlanningHorizon 

Cost(Day) = DelivCost * (delivery(Day) + riskOG(Day, calcday)) 
If Cost(Day) < MinCost Then 

MinCost = Cost(Day) 
End 
Next Day 

'412 Cost rninimum? 

If MinCost <= Cost(l) And MinCost < Cost(PlanningHorizon) Then 

'418 Calculate RAD 

Call CRAD(LHHorizon, SafetyHorizon, RouteHorizon, rnincost, RADR) 
End if 
Next customer 



End Sub 
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Function delivery(FutureDay) 

delivery = CustomerYearlyDemand * (1 / (0.8 * Customer.TankSize) - 1 / (Dem(FutureDay))) 
End Function 

Function riskOG(today, FutureDay) 
'Rho is the use sepecified ratio of cost og OG to cost of delivery 

riskOG - normdist(FutureDay) - normdist( today) / (1 - normdist(today)) * Customer.Rho 
End Function 

Function Dem(FutureDay) 

Dem = mu* (FutureDay - lastdelivery) + hf * (DD(FutureDay) - DD(lastdelivery)) 
End Function 



Figure 2C 
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Figure 3a: Step 120 Risk Estimation 
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Figure 3b: Corresponding pseudo code (When Applicable comments point to the steps in the 
flowchart) 

Sub CRAD (LHHorizon, SafetyHorizon, RouteHorizon, mincost, RADR) 

Dim Threshold As Double 
Dim CostTolerance As Double 

Dim LastCustomer, LastSeq, Customer, Pos, Day, Span, RadSeq As Integer 
Dim RADN(), RADR(RouteHorizon) As integer 
Dim PlaniiingHorizon as Integer 
Dim Day, Dayshift as Integer 
Dim Schedl, Sched2 as Integer 

Dim Db As Database 

Set Db = CurrentDb 

'500 Read Safety horizon 

PlanningHorizon = RouteHorizon+SafetyHorizon+LTHorizon 
Redim RADN(PlanningHorizon) 

Threshold = 0.9 

CostTolerance = MinCost/Threshold 

'502 Calculate RAD 

For Day = 1 To PlanningHorizon 
RADN(Day)=0 

Cost = DelivCost * (delivery(Day) + riskOG(Day, calcday)) 
If Cost < tolerance Then 
If (WeekDay(Now + Day) o 1 And WeekDay(Now + Day) o 6) _ 

Or Now + Day < RouteHorizon _ 
Then 

RADN(Day) = 1 
End If 
End If 

If day > RoutingHorizon+safetyHorizon then 

Priority=2 
Else 

Priority=l 
End If 
Next Day 

'504 Calculate adjusted Min RAD & Max RAD 

For day=l to planninghorizon 
If RADN(day)=l then 
Dayshift=day-saferyhorizon 
If dayshift > routehorizon then 

Day=planninghorizon 
Else 

If dayshift <=0 then 

RADR(1)=1 
Else 

RADR(dayshift)=RADN(day) 
End if 
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End if 
Next day 
Schedl=0 

Forday=l to routehorizon 

Schedl=Schedl+RADR(day) 
Next day 

If Schedl=0 then delete 
Sched2=0 

'510 Adjust for specific day of the week 

For day=l to routehorizon 

IfRADR(day)=l & Customer.FIX(Day)=l then 

RADR(Day)=l 
Else 

RADR(Day)=0 
End if 

Sched2=Sched2+RADR(Day) 
Next Day 

For day=l to routehorizon 

If Customer.FIX(Day)=l & Schedlo 0 & Sched2=0 then 

RADR(Day)=l 
Next Day 
End Sub 



Figure 3 c 
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Figure 4a: Step 130 Routing 
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Figure 4b: Corresponding pseudo code (When Applicable comments point to the steps in the 
flowchart) 

Sub CRMRoute() 

Dim LastCustomer, LastSeq, Customer, Pos, Day, Span, RadSeq As Integer 

Dim RAD As String 

Dim FirstDay, LastDay As Integer 

Dim Demand, AvrDemand, Demand 1, Demand2, CapacityRatio() As Double 
Dim Cust, Count As Integer 

Dim Db As Database 

Dim Customers, RadList As Recordset 

Set Db = CurrentDb 

'600 read information 

Set Customer=Db.openrecordset(Modelreg) 

'602 Select first customer 

Move Customer.first 
Move Customer.last 
CustomerCount=Customer.Recordcount 
Move Customer.first 

For Cust=l to CustomerCount 
RAD = "" 

'604 Construct routing string based on RAD 

For Day = 1 To RouteHorizon 

If RADR(Day)=l then 
RAD = RAD & Day 

End If 
Next Day 

'606 File routing string 

If RAD o Trien 
.Edit 

Customer! RAD = RAD 
.Update 
End If 

Move Customer.next 
Next Customer 

' Customer demand assignement 

FirstDay - Val(Left(RAD, 1)) 
LastDay = Val(Right(RAD, 1)) 
Demand = Dem(FirstDay - 1) 
Count = 0 

For Day — FirstDay To LastDay 
If Mid(RAD, Day, 1) o Then 
Demand = Demand + consumption 
Count = Count + 1 
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End If 

Next Day 

AvDemand = Demand / Count 

'614 Capacity adjustment 

For Day = 1 To RouteHorizon 
Demand 1 = 0 
Demand2 = 0 
Customers.MoveFirst 
For Cust - 1 To LastCustomer 
Span = length(RAD) 
For Pos = 1 To length 
If Val(Mid(Pos, 1) = Day) Then 
Demand 1 = Demandl + Dem(Day) 
Demand2 = Demand2 + AvDemand 
End If 
Next Pos 
Next Customer 

CapacityRatio(Day) = Demand2 / Demandl 
Next Day 

'618 Assign appropriate routing string 

Db.TableDefs.Delete "RADList" 

DoCmd.OpenQuery ("SetUpRADList") 

Set RadList = Db.OpenRecordset("RADList") 

RadListMoveFirst 

RadList.MoveLast 

RadSeq = RadList. RecordCount 

For Seq = 1 To RadSeq 

Span = length(RadListlRAD) 

For Pos = 1 To RadSeq 
Day = Mid(RadList!RAD, Pos, 1) 
truck(Day).specialty = RAD 

'620 Adjust capacity 

truck.(Day).capacity=Truck(Day).capacity*Capacity ratio(day) 
Next Pos 
Next Seq 

'626 Routing 

Call OVRalgorithmO 

End Sub 



Figure 4c 
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Figure 5a: Step 310 Operational audit procedure 
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Figure 5b: Corresponding pseudo code (When Applicable comments point to the steps in the 
flowchart) 

Sub crecon() 

' 700 Read input information 
Start: 

702 Select First customer 

Do Until EOF 

Read "Degree days" 
Date 

DegreeDay 

CumulativeDegreeDay 

Read "Last delivery spec" 
CustomerNumber 
TankS ize 
LastDeliveryDate 
LastDeliveryGallons 
LastPercent 
PreviousDeliveryDate 
PreviousPercent 

Read "Forecasting Parameters" 
CustomerNumber 
DailyUse 
HeatFactor 
Lag 

Residual (Month) 
Seasonal (Month) 
SeasonalProfile (Month) 
DailyStandardDeviation 
CountDeliv 

' Initialization 

GapBound =15 

level = TankSize * LastPercent 

Safety = 0 

Now = LastDeliveryDate 

Gap = InvNorm(l / (CountDeliv * GapBound)) 

704 Calculated planned delivery size 

For Day = PreviousDeliveryDate To LastDeliveryDate 
level = level - SeasonalProfile(Month(Day)) * _ 

(DailyUse + HeatFactor * DegreeDay(Day) + Lag 12 * Residual(Month(i))) 
Next Day 

Safety = Gap * DailyStandardDeviation * __ 

sqt(Sum(i = PreviousDeliveryDate, LastDeliveryDate)(Seasonal(Month(i)))) 

706 Extract Basic delivery statistics 
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'Adjust actual delivery 

LastDeliveryGallons = LastDeliveryGallons * 0.8 / LastPercent 



'Calculate days to runout 
DaysToRunOut = 0 
StatisticalOG = level 
Day = LastDeliveryDate 
Do Until StatisticalOG < 0 
DaysToRunOut = DaysToRunOut + 1 

StatisticalOG = StatisticalOG - SeasonalProfile(Month(Day)) * _ 

(DailyUse + HeatFactor * DegreeDay(Day) + Lag 12 * Residual(Month(Day))) - 
Gap * DailyStandardDeviation * _ 

sqt(Sum(i = PreviousDeliveryDate, LastDeliveryDate + DaysToRunOut) _ 
(Seasonal(Month(i)))) 
Loop 

* 708 Store performance information 

write 
CustomerNumber 
DaysToRunOut 
level 
Safety 

LastDeliveryGallons 
710 Last Customer? If no then 712 move to next customer 
Loop 

'if yes then stop 
End Sub 



Figure 5c 
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Figure 6: Step 320 Overall flow of the modeling step 
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Figure 7b: Corresponding pseudo code (When Applicable comments point to the steps in 
flowchart) 



Sub clusteringO 

' 900 Select all statistically predictable customers 

Read "CompletedProcessing" 
CustomerNumber 
TankSize 
RoutingMethod 
ConsumptionType 
Seasonal profile 
StartDate 
CountDeliv 
TotalGallons 
MeanGallons 
StdGallons 
Meanlnterval 
Stdlnterval 
DailyUse 
HeatFactor 
Lagl2 

Seasonal (Month) 
Residual (Month) 
DailyStandardDeviation 
YearConsumption 

' 902 Cluster tanks for each consumption type 

Define TankThreshold() 

AverageConsumption = DailyUse + HeatFactor * YearDD / 365 
Sort by ConsumptionType tanksize output SortedTanks 

For Each TankSize In ConsumptionType 

If Not First.TankSize Then 
Gap = TankSize - Lag(TankSize) 

End If 
Next TankSize 

sort by ConsumptionType and descending Gap 

Do Until Gap < TankThreshold 

TankBound = TankSize 

output ClusterBoundaries 
Loop 

f 904 estimate distribution of consumption for each cluster 

Cluster = First(ClusterBoundaries) 

For Each TankSize In ConsumptionType 
If TankSize < TankBound Then 
Percentile (Consumption) 
Else 

CDailyUse = SafetyLevel(DailyUse) 
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CHeatFactor = SafetyLevel(HeatFactor) 
MoveNext ClusterBoundaries 
End If 
Loop 

'906 estimate average consumption for each consumption type 

for each ConmsumptionType 

MDailyUse = mean(DailyUse) 

MHeatFactor = mean(HeatFactor) 
Loop 

' 908 Sort coefficient of variation for each consumption type 

For Each TankSize In ConsumptionType 

CV = StdDev / Consumption 
sort by consumptiontype consumption CV 

End Sub 



Figure 7c 
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Figure 8a: Step 808 Matching step 
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Figure 8b: Corresponding pseudo code (When Applicable comments point to the steps in the 
flowchart) 

Sub MatchingO 
Do Until EOF 

* 1000 select customer without model 

Read "PostponedProcessing" 
CustomerNumber 
TankSize 
RoutingMethod 
ConsumptionType 
SeasonalProfile (Month) 
StartDate 
CountDeliv 
TotalGallons 
MeanGallons 
StdGallons 
Meanlnterval 
Stdlnterval 
DailyUse 
HeatFactor 
Seasonal (Month) 
Residual (Month) 
YearConsumption 

'Clustering criterion for variance 

Define VarboundO 

' 1002 Tank size known 

If TankSize - 0 Or TankSize = Null Then 

' 1004 Use average consumption parameters 

Select Case ConsumptionType 
Case "Heat" 

DailyUse = 0 

HeatFactor MHeatFactor 
Case "Non-Heat" 

DailyUse = MDailyUse 

HeatFactor = 0 
Else 

DailyUse = MDailyUse 
HeatFactor = MHeatFactor 
End Select 

'1006 Match tank size to nearest cluster, associate corresponding _ 
consumption parameters 

Else 

Select Case ConsumptionType 
Case "Heat" 
For Each TankBound In ConsumptionType 
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If TankSize < TankBound And TankSize >= Lag(TankBound) Then 

DailyUse = 0 

HeatFactor = CHeatFactor 
End If 

Next TankBound 
Case "Non-Heat" 
For Each TankBound In ConsumptionType 
If TankSize < TankBound And TankSize >= Lag(TankBound) Then 
DailyUse = CDailyUse 
HeatFactor = 0 
End If 
Next TankBound 
Else 

For Each TankBound In ConsumptionType 
If TankSize < TankBound And TankSize >= Lag(TankBound) Then 
DailyUse = CDailyUse 
HeatFactor = CHeatFactor 
End If 
Next TankBound 
End Select 
End If 

1 1008 estimate seasonal characteristics 
Lagl2 = 1 

For Month = 1 To 12 
If Seasonal o Null Then 

Residual = Seasonal / 30.4 * (DailyUse + HeatFactor * YearDD / 365) 
Else 

Residual = 1 / 30.4 * (DailyUse + HeatFactor * YearDD / 365) 
End If 
Next Month 

'Variance estimation 

EstConsumption = DailyUse + HeatFactor * YearDD / 365 
Read "VarCluster" 

'1010 Cluster customer nearest estimated consumption parameters 

For Each Consumption In VarCluster(ConsumptionType) 
If Consumption <= EstConsumption + Varbound And _ 
Consumption > EstConsumption - Varbound Then 
Percentile (CV) 
End If 
Next Consumption 

'1012 estimate distribution of coefficient of variation of the cluster 
estCV - Safety(CV) 
'1014 estimate variance 
DailyStandardDeviation = CV * Consumption 



Figure 8c 



Docket No.: RMR-0003 

App No.: Not yet assigned Filed: Herewith 

* Title: Method for Forcasting Consumption and Generating Optimal Delivery Schedules for 

Vehicles Involved In Delivering Propane and Other Consumables to End Consumers 
Inventors: Dominique M. Richard and Joel Adler 

Attorney: Michael P. Dunnam Phone: (2 1 5) 568-3 1 00 

Sheet 21 of 21 



write "CompletedProcessing" 
RC 

CustomerNumber 

TankSize 

RoutingMethod 

ConsumptionType 

SeasonalProfile (Month) 

StartDate 

CountDeliv 

TotalGallons 

MeanGallons 

StdGallons 

Meanlnterval 

Stdlnterval 

DailyUse 

HeatFactor 

Lag 12 

Seasonal (Month) 
Residual (Month) 
DailyStandardDeviation 
YearConsumption 

'1018 move to next customer 

MoveNext 
'1016 last customer 
Loop 
End Sub 



Figure 8d 



